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' Abstract 

The time dependent Schrodinger equation inclusive of curvature effects is developed for a spinless 
electron constrained to motion on a toroidal surface and subjected to circularly polarized and 

Qh. linearly polarized waves in the microwave regime. A basis set expansion is used to determine 

c/3 , 

c/2 \ the character of the surface currents as the system is driven at a particular resonance frequency. 

C^ ' 

O ' Surface current densities and magnetic moments corresponding to those currents are calculated. It 

CJ ' is shown that the currents can yield magnetic moments large not only along the toroidal symmetry 

C/2 . 

J>^, axis, but along directions tangential and normal to the toroidal surface as well. 
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1. INTRODUCTION 

The control of a nanostructure's state through electromaenetic interactions is of funda- 
mental and practical interest [103113. Considerable effort has been d.reeted towards the 
study of flat quantum rings, which, because of their topology, giv e rise to Aharanov-Bohm 
and persistent current effects 0, fl 3 B Q, H Q, H H H H H ^ 



Currents 



generated on a ring are restricted to yield magnetic moments perpendicular to the plane of 
the ring; however, it is conceivable that structures with different topologies are capable of 
producing magnetic moments in different directions, a feature potentially useful when such 
structures are employed as nano-device elements. In this work, a toroidal structure in the 
presence of an electromagnetic wave on the microwave scale is considered. In contrast to 
rings, it is possible to excite modes for which currents circulate around the minor radius of 
the toroidal surface T^ (called in what follows "solenoidal modes" since the focus will be on 
magnetic moments) in addition to the modes in which currents circulate azimuthally around 
the major radius ("dipole modes"). 

The approach adopted here will be to develop 

H = Ho + V{t) (1) 

with Hq the Hamiltonian for a spinless electron constrained to motion on T^ inclusive of 
surface curvature effects 2l|, and V{t) the time dependent electromagnetic interaction with 
the form appropriate to circularly and linearly polarized microwaves propagating along the 
symmetry axis of the torus (here the z-axis) at a particular resonant frequency of the system. 
It has been shown that for the case of no external vector potential, it is possible to treat 



the normal degree of freedom as a spectator variable [2^. However, inclusion of a vector 
potential A with components normal to the surface complicates the situation. Reducing 
the full three-dimensional problem to a two-dimensional effective model serves to simplify 
matters considerably. The resulting Schrodinger equation is then solved with standard 
methods by taking [h = 1) 

^(t) = J2(^nuit)Xn{0)exp[iiy(l)]exp[-iEnut] (2) 

ni/ 

with the Xn{d)exp[—ih'(j)]exp[—iEnut] eigenstates oi Hq. Surface current densities (SCDs) are 
determined and magnetic moments are calculated, from which it is shown that is it possible to 
selectively generate dominant oscillatory magnetic moments in different directions depending 
on the character of the applied signal. 

The remainder of this paper is organized as follows: in section 2 the Hamiltonian on T^ is 
developed by beginning with a three dimensional formulation and proceeding to restrict the 
particle (accounting for curvature)to T^. In section 3 the eigenstates of Hq that comprise 



the basis set are given and the solution method detailed. Section 4 presents results in the 
form of plots of time dependent SCDs and magnetic moments. Section 5 is reserved for 
conclusions. 



THE SYSTEM HAMILTONIAN 

This section presents a derivation of H{t) inclusive of surface curvature (SC) effects. The 



necessity of accounting for SC has been shown elsewhere 2^. The central idea is upon 
restricting a particle in the neighborhood of a surface to the surface proper, a geometric 
potential Vc is necessary to capture the full three dimensional spectra and wave functions 
with good fidelity. 

In this work, Vc will not be the only geometric potential. The normal component of A 
also couples to SC |l6|, |23|. It should also be noted that in general, the Laplacian on the 
surface must be modified due to curvature. In the case of T^ the symmetry of the system 



precludes such modifications [21 1. 



The Schrodinger equation for a spinless electron in the presence of a vector potential A(t) 
(reinserting h) is 

H=^ y + eAit))^ = ^h^-. (3) 



2m \i J dt 

To derive the gradient appearing in Eq. (3), let 6^,6^,62 be unit vectors in a cyhndrical 
coordinate system. Points in the neighborhood of a toroidal surface of major radius R and 
minor radius a may then be parameterized by 

x(6', 0, g) = (i? + a cos6')ep + a sm.9ez + ge„ (4) 

with e„(to be given momentarily) everywhere normal to T^, and q the coordinate measuring 
the distance from T^ along e„. The differential line element rfx is then 

(ix = a(l + kiq)eed9 + W{9, g)(l + k2q)e^d(f) + e„(ig (5) 

with 

e„ = cos6'ep + sin6'ez (6) 

eg = —siiaOep + cos6'ez (7) 

W{e,q) = R+{a + q)cose (8) 

and the toroidal principle curvatures 

h = - (9) 

a 



The result for V is then 

1 d 1 d d_ 

^ ~''\a + q)de'' ''V(^, q) 90 "^ ^"^ dq ^ ' 

To proceed with the development of H{t) specific choices must be made for A(t). Here 
the vector potential corresponding to a CPW 

A(t) = — [ — sui{kz — uotyi + cos(A;2; — uot)'^ (12) 

UJ 

with Eq the field amplitude, uj the frequency and z = [a + g)sin6' is taken. The LPW case 
can be obtained trivially from Eq. (12). 

The procedure by which Vnis obtained from Eos. f3)-(ll) is well known (more detailed 



derivations can be found in [2J, |2^ |2a, \2J\, |2^ 1221, 120, lM|) and is summarized briefly below. 



Begin by making a product ansatz for the surface and normal parts of the wave function 

{l + 2qh{e) + q^k{e))2 
with h{6), k{9) the mean and Gaussian curvatures of the surface |3J|, and imposing conser- 
vation of the norm through 

\m{e,cl>,q)\\l + 2qh{e) + q^k{e))<n:dq^\Xs{0.<t>)?\XN{q)?(n:dq (14) 

with (is the differential surface area. The Laplacian appearing in Eq. (3) takes the form 
(the subscript on the gradient operator indicates only surface terms are involved) 

V=^V| + 2M^)| + ^ (15) 

with q appearing in V| as a parameter that can be immediately set to zero. The differenti- 
ations in Eq. (15) act upon the wave function given in Eq. (13) to yield 

in the g ^ limit. The /i^(6') — k{9) term above is proportional to Vc{0)., and along with 
V/v(5'), the remaining q differentiations can produce solutions confined arbitrarily close to a 
thin layer near the surface S. 

As noted earlier, in addition to Vc there is a second geometric potential that arises 
from the normal part of the A ■ V term in Eq. (3). Because any effective potential on S 



should only involve surface variables the operator AN{6,(j),q)^ must be addressed. One 
prescription to deal with this operator has been given in jl6| and another in 23| which is 
chosen here. Begin by noting that the energy level spacing corresponding to the eigenvalues 
for the Xn{(i) normal term is large compared to surface eigenvalues so it is safe to assume 
there is negligible mixing among the Xn{(i) (this assumption was shown to be a reasonable 
one in J2^). Let G = 1 + 2qh{6) + q^k{6) and proceed to integrate out any g-dependence by 
writing 



XN{,q) 



Ajv(^,0,g) 



d XNJq) 
dq G^l^ 



Gdq 



(17) 



'0 GV2 

Given any well behaved A^y it is simple to establish after an integration by parts that 
the resulting effective potential V"^°'^ is (with constants appended) 



V^'i 



ieh 
m 



h{9)AN{9,(P,0) + 



ldAM{9,(P,q) 



dq 



q=0. 



(18) 



The ensuing Schrodinger equation can put into dimensionless form by defining 



a = a/R 
F{9) = l + a cos9 

In = — 
e 

after which Hq may be written 

(9^ a sin 9 d c? d"^ 



Hn 



+ 



d^9 F{9) 09 F^ (9)8^(1) 
The interaction potential V{t) with V"^""^ included becomes 



+ h^{9)-k{9). 



V{t) 



sin^(sinficos0 — cosr2sin0)-7— + 

a 09 



2ia'^aTEQ 

1 Id 

-— — (sinOsin0 + cosf] cos0) — — 



1 d 



where 



1 + 2a;cos6' ~ 
^2^FW 

ANiOA) = 



— sin2^(cosi7cos0 + sinf2sin0) 



-sinr2cos6' cos + cos r2cos6' sin 0, 



(19) 



(20) 
(21) 



VL{9) = ka sm9 — oot 
and the |A(t)p term has been neglected. 



(22) 



Note that the length gauge jl| is not employed here. Although implementing the gauge 
transformation in a curved geometry that leads to the dipole term does not involve any 
inherent difficulty, it is useful to report expressions for matrix elements should an extension 
to this problem involving ionization be desired (as would be possible for finite layer rather 
than ideal surface confinement), or for cases involving thin tori with minor radii on the order 
of a few angstroms subjected to visible laser light. 

3. BASIS SET AND METHOD 

In this work R will be set to 500A, a value in accordance with fabricated structures 



33|, |3J, |35|, |36| , and a = .5, a value which serves as a compromise between smaller a where 
the solutions tend towards simple trigonometric functions and larger a which are less likely 
to be physically realistic. 

The eigenstates of Hq are found by diagonalizing Hq with a 60 state basis set expansion. 
The sixty states comprise five azimuthal functions {u = —2, ..., 2) multiplied by six positive 
and six negative 6 parity Gram-Schmidt (GS) functions. The GS functions are constructed to 



be orthogonal over the integration measure F{6) = 1 + a cos^ 37|. The energetically lowest 
six states are dominated by the constant or cos^ mode. To create a non-zero net current 
around the minor axis of the torus, some part of the wave function which behaves as e*"^ 
must appear; it is not until the seventh state (as ordered by energy) that a sine term appears 
which motivated the choice of the signal frequency as un = u-i—uji. The resonance frequency 
uJij (in natural units) between any two levels i and j is numerically 3.09 x 10^*^ {j3j — /5j)A~^. 
The associated resonance wavelength between the ffist and seventh states corresponds to 
approximately 2 cm. The electric field value chosen here is approximately 10 V/m, a value 
large enough to induce effects but small enough to justify the neglect of the quadratic vector 
potential term in Eq. (3). The signal was applied to the structure for ten periods of the 
inverse resonance frequency, for a time of 2.04 x 10^°r units, which is equivalent to tf = .68ns 
with r = 3 X lO^^ns. 

The matrix elements (nz/|V(t)|?2z/) can be evaluated analytically in terms of Bessel func- 
tions J38|. The system of seven ffist-ordered coupled equations is then solved for the dnu{t)- 
Selection rules between the v and v are sufficient to render the system of coupled equations 
relatively sparse allowing the system to be solved by standard methods. 



4. RESULTS 

In this section SCDs as computed from 

J(t) = — /m[^*(^,0,t)V5^(^,0,t)] (23) 

lit 

are found. Again, the subscript on the gradient operator indicates only surface terms are 
involved. From the J(t), magnetic moments 

1 r2w /•27r 



M(t) = - / / r X 3dA (24) 

2 Jo Jo 

with dA = aW{6)d6d(j) are presented. It is possible to write M(t) in terms of unit vec- 
tors e0,e^,en, but because the interest lies in comparing solenoidal to dipole modes, the 
components along e^, e^, e^ are shown instead. Explicit forms for the magnetic moments 
Mp{t), M^{t), Mz{t) are given in the Appendix. 

Before presenting results for the currents and moments, it is worth showing representative 
plots of the time dependent coefficients di{t). First let L serve as a collective index for the 
values i = 1....6. In Fig. 1, |(ii(t)p and |(i7(t)p for the LPW are shown, and Fig. 2 gives 
the same for the CPW. In both situations, |(i7(t)p is small as compared to the remaining 
|(ii(t)p, but oscillate at a much slower frequency than the |(ij^(t)p. As was noted previously, 
it is di(t) that must multiply the sine terms to combine with the positive 6 parity parts of 
the total wave function to yield currents around the minor radius, so it is the time scale of 
di{t) that will set the time scale and magnitudes of J(t) and M(t). 

Results for SCDs as functions of 6 evaluated at t/ are shown in Figs. 3 and 4 for the 
LPW. Fig. 3 plots Jgitf, 6), the current calculated from employing only the diit), and Fig. 
4 plots Jg{tf,6), the current resulting from inclusion of di{t) summed over the d^lt), both 
at = 0. The results illustrate that the net current resulting from Jg(tf,9) is zero, but 
is non-zero for Jg{tf,6). Figs. 5 and 6 plot the same quantities for the CPW. Results for 
the azimuthal current J,j>{tf, 0) at ^ = and 6 = n are shown in Figs. 7 for the LPW and 
in Fig. 8 for the CPW. The contribution from terms proportional to dj(t) is always much 
smaller than the contribution arising solely from the dilt) terms, so they have not been 
shown separately. No net azimuthal current results from the LPW, and only a small current 
is present in for the CPW case. 

In Fig. 9 the Mj(t) for the LPW are shown for a duration of tj at = 7r/2, ^ = 0. The 
surprising result here is that M^{t) is generally an order of magnitude larger than the dipole 
moment Mz{t). Although not shown here, it is observed that as the torus is traversed in 
the direction, the direction of the largest component of the magnetic moment is consistent 
with a magnetic moment perpendicular to the direction of polarization. The corresponding 



results for the CPW presented in Fig. 10 indicate that the dipole mode here is four orders 
of magnitude larger than any solenoidal modes, a result that holds true at every angle. 

The above results were applicable to an ideal situation wherein the structure is situated 
on a surface transparent to microwaves. It is interesting to consider what is perhaps the 
more realistic case of the incident wave interfering with its reflected part. 

In order to see whether the interference of two polarized waves propagating along the 
z-axis can enhance transitions to solenoidal modes, an incoming electromagnetic wave being 
reflected at a completely reflecting flat surface ('mirror') parallel to the horizontal symmetry 
plane of the torus was considered. The mirror is positioned at variable distances zq from 
the central plane in order to study how superposition of the incoming and reflected linear 
polarized wave fronts might affect the relative amplitude of the solenoidal modes. The 
general arguments will be made for the CPW; the LPW results are a special case of those 
for the CPW. 

At the mirror, the electric field must satisfy the boundary condition of no total transverse 
modes 

EUzo,t) + Eout{zo,t)=0, (25) 

The location and effect of the reflecting mirror on the incoming wave front in this context 
may be simply modeled in form of an additional phase 0o = —^kzo in the reflected field 
amplitude. The signs of uj and zq can be chosen according from the incident direction of the 
incoming wave. From Eq. (25), the reflected wave amplitude must be: 

Eiout{z, t) = Eo[ — cos{kz + Ljt')i + sm{kz + a;t')j] (26) 

with t' = t — to; and to = —2kzQ/uj. Thus, the reflected wave amplitude can be obtained 
from the incoming wave amplitude by letting 

u -^ -u, t-^t', j ^ -j. (27) 

in Eq. (26). 

Several different positions for the position of the mirror plane {zq = a, |, | and |), were 
investigated. The general outcome is simple oscillatory behavior with different but small 
amplitudes while \di\ stays close to 1. For zq = j, all coefficients ^2, . . . , (i? are numerically 
nearly equal to zero, while di stays constant and identical to 1. This is due to the fact 
that the electric field amplitude E reaches its maximum value, thus the field gradient is 
essentially zero over the whole torus. In all of these cases, the system stays essentially in 
the ground state and no transitions are observed. 

I 1 2 

The optimal mirror position that maximized the value |(i7(t)| was determined to be at 

zo ^ 0.0027A, (28) 



where it was found that \dY\ reaches its maximum amphtude of roughly 7 * 10"'''. This is 
an increase of up to 3 to 4 orders of magnitude for positions around j discussed above. 

The analogous results for the interference cases to those presented above are given below. 
Figs. 11 and 12 show |(ii(t)p and 1^7 (t)p plots for the LPW and CPW. It is worth noting 
that while the |(i7(t)p time scale does not vary greatly due to interference as compared to 
the results in Figs. 1 and 2, the |(ii(t)p time scale is very different in both cases. 

The SCD plots for the LPW do not appear qualitatively different than the results already 
presented in Figs. 3, 4 and 7, so in the interest of conciseness those results are not presented 
here. However, an interesting result is shown in Fig. 13, wherein a clear signature of a 
circulatory Jg(tf,9) SCD at = is shown for the CPW due to the presence of the dT{t) 
term. Fig. 14 shows there is also a circulatory azimuthal SCD at ^ = tt. 

The LPW magnetic moment results are given in Fig. 15. The frequency of oscillation 
has again been reduced, and the overall magnitude of Mz{t) is smaller than that of M^{t). 
Fig. 16 shows that for the CPW, Mz{t) > M^{t); this is observed to be independent of the 
azimuthal point at which the Mj(t) are calculated. 

5. CONCLUSIONS 

This work presents a general framework for calculating surface current densities and mag- 
netic moments on a toroidal surface in the presence of an electromagnetic wave. A proof 
of principle calculation demonstrating that polarized microwaves can cause circulating sur- 
face currents around the toroidal minor radius in addition to azimuthal currents was given. 
Rather than employ the dipole approximation, matrix elements of the electromagnetic term 
were evaluated in closed form. While there can be little doubt that the dipole approximation 
would be perfectly valid for the system considered, there are cases involving ionization or 
surface states with energies characterized by much smaller minor radii that warrant exact 
expressions. 

The magnetic moments calculated in this work show that it is in principle possible to 
manipulate surface currents in a manner that causes the moments to "point" predominantly 
in particular directions at certain times. Our preliminary results indicate that interference 
due to reflection could potentially play an important role in the development of this topic. 

The realization of the model problem considered here may be physically realizable with 
current experimental methods by placing a thin layer of a good conductor over an InAs 
toroidal structure. The extension of this work to metallic carbon nanotube tori is likely 
possible but would require some effort to account for the lattice. Preliminary work indicates 
that modelling the carbon sites by weak delta function potentials requires a larger basis set 



than employed here. 
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Appendix 

This appendix presents closed-form expressions for the integrated magnetic moment and 
its components. The magnetic moment M(t) has been deferred as: 



M(t) 



27r /•27r 



^0 



r X J{r, t)dA 



with the radius vector r on the T^ surface generated by 

r = {R + a cos 6)ep + a sin Oe^ 

and the surface gradient being 

^ I d Id 

V5 = 60 - 7- + e, 



a 86 R + a cos 9 (90 



(A.l) 



(A.2) 



(A.3) 



The time dependent current density J(r, t) in Eq. (A.l) is defined in Eq. (23) with the wave 
function ^(r, t) given in Eq. (2) as solution of the time dependent Schrodinger equation. 
Here express the 9 dependent part of the eigenstates, Xn{9), directly in terms of sums of 
cosines with coefficients C^^ for positive 9 parity solutions, or as sums of sines with coefficients 
D'^ for negative 9 parity solutions respectively: 



Xni9) 



-^Y.CnCos{k9) for n = 1...6; 



fc=0 

5 



(A.4) 



-^T.Di sm{k9) for n = 7. 



k=l 



Now introduce the following constants Mq, C and A: 

Mo = ^R\ C '^ 



Arrip 



2mp 



A = 7rR\ 



(A.5) 



Applying Eq. (A.2) through Eq. (A.5) the time dependent magnetic moment M(t) in Eq. 
(A.l) can be written as integral via the surface variables 9 and with F{9) = 1 + acos^: 



2-77 TV 



M(t) 



-iTrMn 



d9 



aF{9)J^e^ - F (9) (cos 9 + a) Jge^ - a'^sm9J^ep (A. 6) 



with following expressions for the current densities: 



9^ 9^* 

J^(r,t) = M/*(r,t) — (r,t)-vl/(r,t)^(r,t); 



(A.7) 
(A.8) 
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The currents in Eq. (A. 7) and Eq. (A. 8) can be re-expressed as 



(A.9) 
(A.IO) 



*j 



with the time independent current expressions J^J {6, 0) and J^-^ {9, 0) presented in the fol- 
lowing way (cjjj = -jf- 



Ei — Ej 



Uii = Ui - Ua : 



^3 — h — h 1 ""ij — ""i '^3/ 



2Q-^''i3'f>ly*^^ - -^* 



06 ^^'oe 



(A.ll) 
(A.12) 



For the final results for the magnetic moments we prepare the integrations over 6 and 0. 
The 0-integration yields following three integrals: 



2n 



)e, = 27i6{ui - z/,-)e'"'^*e„ 




27r 



2tt 



,^i{u>^jt-u,,4>)Q^ = TT 



^e'^'''^^-'''^'t''^e, = TT 



Sp + Sn)ey + i[6p - 6n)e^ 



Op -\- On I Gj. '■[Op On I Gy 





(A.13) 




(A.14) 


e*^'^-*, 


(A.15) 



which includes the Kronecker delta expressions 



6p = 6{uj -i^i + 1); Sn = S{uj -Vi- 1); 

5^ = 5{k-l); 6+ = 6{k + l); 

6^^ = 6{k-l± 1); 6^^ = 6{k + l± 1); 

6^^ = 6{k-l±2); 6^^ = 6{k + l±2). 



(A.16) 
(A.17) 
(A.18) 
(A.19) 



The final results for the ^-integration are based on following integrals of products of sines 
and cosines: 



a 



deaF{e)cos{ke)cos{W) = - 



d9aF{9)sm{k9)sm{l9) 



a 



6^ + 6^ + ^(61' + 6+' + 6-' + 6+' 



6_-6+ + ^(6z' + 6l:'-6-'-6X' 



(A.20) 
(A.21) 



13 



d9 F{9) (cos 9 + a) cos{k9) cos{l9) 



^a ((5_ + 5+) + ^ (i + a') (r 1 + 61' + 5;^ + 5+^) + ^ ( r ^ + 5+2 + 51' + 5^' 



(A.22) 



rf^ F{9) (cos ^ + a) sin(A;^) sin(/^) 



:^a 6^-6,]+;- l+a'] 6z' + 6l'-6-'- 6l' ] + :^ 5z' + 61' - 51' - 5 



a 



:-2 , r+2 r-2 r+2 



/2 
rf^ a^ gjj^(^) cos(A;^) sin(/^) = ^ (5-' - 5+^ - r ^ + 6 



(A.23) 
(A.24) 



These results in Eq. (A. 20) to (A.24) combine to the final ^-integrals of products of eigen- 

states Xi, Xf 



k,l=0 



(5_ + (5+ + f (r 1 + 61' + 6-' + 5+1) 



d9x^{9)xA0)Fi()) = ^{ 



5 



E D'D^, 



^l x-^ I x+^ x-^ A+i 



(5_-5+ + f(r^ + r^-5;^-5 



for i, j = 1 . . .6; 

for i = j = 7; 

(A.25) 



rf^(x.^-X,^)i^W(cos^ + a) 



5 5 



A;=0 /=1 



^EE^- r{5.+6,) + l{l+a'){6-J + 6l' + 5-' + 6X' 



+1 



+^{sz' + Sl' + 6-' + 6X' 
^a (5_ -S+)+l{l + a')[ 6-J + 5+^ - S-' - S^ 



a 



+:_{Sz' + 6l'-6l'~6+' 



C^Dj for i = 1 . . .6 and j = 7; 

{-C^D^i) ior i = 7 and j = 1 ... 6; 



(A.26) 
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ddUd)xAd)sind = 

1 v^ v^ / 1 ^1 1 ^^\ { CfD\ for i = 1 ... 6 and J = 7; 

4^ t^ ^ ^ ^1^ C^D\ for i = 7 and J = 1 ... 6. 

(A.27) 

In the final expression for the magnetic moment M, re-express e^ and e^, in terms of e^ and 
e^ while introducing the fixed angular variable 0: 

e^ = — sin0e<^ + cos0ep; 

e^ = cos0e<^ + sin0ep. (A. 28) 

The ^-integration finally yields three contributions Mi(t), M2(t) and M3(t) for the total 
magnetic moment. The dipole mode Mi (= M^) parallel to the central symmetry axis can 
be determined using the definitions (if (t) = Re {di{t)) and (if (t) = Im {di{t)) as: 

27r 7r 

Mi(t) = M,{t) = -ianMo f d^ f d9F{9)J^{9,(P,t) e, = 

-TT 

7 

= e^ (4a7rMo) ^ (z/j + Uj) S {vi - Uj) 

df{t)df{t) + d{{t)d]{t)\ cos{ujijt) - Uf{t)d]{t) - dl{t)df{t)\ sm{ojijt) 

dexiiO)XiiO)F(ey, AaTTMo = -^{2TraR). (A.29) 

2me 

The two solenoidal modes are expressed by M2(t) and M3(t). M2 is proportional to the 
integrated current term Jg in Eq. (A. 8): 
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27r TV 



M2(t) = iTiMo 



d9F{9){cose + a)Jg{9,(l),t) e^ 



-2nMo Y^ 



dfit)d\{t) - dl{t)df{t)\ (6pCos{uJijt - 0) + 6nCos{uJijt + 

+ [d^{t)df{t) + dl{t)d]{t)^ {dp sin(cj,jt - 0) + 5„ sin(cuijt + q 
+ {df{t)d]{t) - dl{t)df{t)^ {Sp sm{uijt - 0) - 5„ sin(cuijt + 

+ {df{t)df{t) + d{{t)d]{t)\ Up cos{uijt - (f)) - 5n cos{uJijt + 0)) 



^Xj -^Xi 



9^ '"^ a^ 



eh 



dO(x^^ - Xj^)F{e) {cose + a); 2'kM, = ^{-kR^). (A.30) 



2mp 



The second solenoidal mode M3 is proportional to the integrated current J^ in Eq. (A. 7): 

2-77 TT 



M3(t) = iaVMo I d(j) I desm9J^{e,(f),t)e^ 



-W-kMo ^ 



-TT 



[z/i + Z/j, 



*J = 1 



df (t)cij- (t) - rfi it)d^{t) 6p COs{LJijt ~(p)-6n COsiyUJijt + 



5n sin(cjjjt + 



+ df{t)df{t) + 0?,' (t)o?j (t) 5p sin{ujijt 



+ 



- {df{t)d]{t) - dl{t)df{t)^ {dp sin{ojijt - 0) + 5„ sin(cjijt + 0)) 



+ df{t)df{t) + d\ {t)d) (t) 6p cos{uJijt - 0) + 5„ cos{uJijt + 



rf^X.Wx.Wsin^; 



efi 



2aVMo = - — (tto^) . (A.31) 



2771, 



The two components mentioned in the text, Mp{t) and M^{t), can then be very easily ob- 
tained from Eq. (A.30) and (A.31) by adding M2 and M3 and projecting into the directions 



Bp or 60. 
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Figure 1 {Eticitidea aud Jack) 

FIG 1. Ground state |(ii(t)p and seventh state |(i7(t)p probability amplitudes for the system 
subjected to a LPW along the toroidal symmetry axis (time scale here: r = 3 x 10~^^ ns). 
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Figure 1 (Ebciuflfia acbd Jack) 

FIG 2. Ground state |(ii(t)p and seventh state |(i7(t)p probability amplitudes for the system 
subjected to a CPW along the toroidal symmetry axis. 



20 



/ 



0.5- 



-0.5- 



-1 -0.5 0.5 1 

Figure 3 (Encinosa and Jack) 

FIG 3. JgltfjO) in units of Jq = eh/m (longest arrow = .00015) plotted on the torus (in units of 
R = 1) at (/) = in 9 = 7r/4 intervals for the LPW case. The net current integrated over the loop 
is zero (t/ = 0.68ns). 
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Figure 4 (Encinosa and Jack) 



FIG 4. Jg{tf,9) plotted on the torus (in units of R = 1) at (j) = in 9 = it /A intervals for the 
LPW case. The longest arrow is 3 x IO^^Jq. The net current integrated over the loop is non-zero 
and arises from the presence of the negative 9 parity d'j{t) amplitude. 
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Figure 5 (Encinosa and Jack) 



FIG 5. Jg{tf,6) plotted on the torus (in units of i? = 1) at = in = 7r/4 intervals for the 
CPW case. The longest arrow is .034 x 10"^ Jq. The net current integrated over the loop is zero. 
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Figure 6 (Encinosa and Jack) 

FIG 6. Jg{tf,9) plotted on the torus (in units of R = 1) at (j) = in 9 = it /A intervals for the 
CPW case. The longest arrow is 6.2 x IO^^Jq. The net current integrated over the loop is small, 
but non-zero. 
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Figure 7 (Encinosa and Jack) 

FIG 7. Jff){tf,(j)) plotted on the torus (in units oi R = 1) at 9 = in (/) = tt/2 intervals for the 
LPW case. The longest arrow is .OIUq. The net current integrated over the loop is zero. 
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Figure 8 (Encinosa and Jack) 



FIG 8. J(j){t f , (f>) plotted on the torus (in units of i? = 1) at = in </> = 7r/2 intervals for the 
CPW case. The net current integrated over the loop is non-zero. The longest arrow is .IUq. 
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Figure 9 (Encinasa and Jack} 

FIG 9. Magnetic moments Mz{t) (top) and M^{t) (lower) in units of Mq = ehR'^/Am at 4> = 7r/2 
for the LPW case. The longest arrow is .OOQJq- 
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Figure 10 (Eacinesa and Jack) 



FIG 10. Magnetic moments Mz{t) (top) and M^{t) (lower) in units of Mq = ehR'^/4m at 
for the CPW case. 
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Figure 11 {Bnciuosa and Jack} 

FIG 11. Ground state |(ii(t)p and seventh state |d7(t)p time dependent probability amplitudes 
for the system subjected to a LPW allowed to interfere with its reflection off an ideal mirror. 
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FIG 12. Ground state |(ii(t)p and seventh state |d7(t)p time dependent probability amplitudes 
for the system subjected to a CPW allowed to interfere with its reflection off an ideal mirror. 
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Figure 13 (Encinosa and Jack) 



FIG 13. Jgiifj G) plotted on the torus (in units of i? = 1) at = in = 7r/4 intervals for the 
CPW case where the wave is allowed to interfere with its reflection. A circulating SCD results. 
The longest arrow is 1.2 x IO^^Jq. 
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Figure 14 (Encinosa and Jack) 

FIG 14. j\{tf,(j)) plotted on the torus (in units of i? = 1) at ^ = vr in = 7r/4 intervals for 
the CPW case where the wave is allowed to interfere with its reflection, showing that a circulating 
SCD results in the azimuthal direction as well as around the minor radius as per Fig. 13. The 
longest arrow is 1.5 x IO^^Jq. 
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Figure 15 {EncinoaaaiudJack) 



FIG 15. Magnetic moments Mz{t) (top) and M^{t) (lower) in units of Mq = ehB? /Am at 
for the LPW case where the wave is ahowed to interfere with its reflection. 
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FIG 16. Magnetic moments Mz{t) (top) and M^{t) (lower) in units of Mq = ehB? /^m at 
for the CPW case where the wave is allowed to interfere with its reflection. 
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